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Plasmoid structures in fast reconnection in low-beta plasmas are investigated by two- 

^ dimensional magnetohydrodynamic simulations. A high-resolution shock-capturing 

' (-h ' code enables us to explore a variety of shock structures: vertical slow shocks behind 

^ the plasmoid, another slow shocks in the outer-region, and the shock-reflection in the 

^ front side. The Kelvin-Helmholtz-like turbulence is also found inside the plasmoid. 
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It is concluded that these shocks are rigorous features in reconnection in low-beta 
plasmas, where the reconnection jet speed or the upstream Alfven speed exceeds the 
sound speed. 
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I. INTRODUCTION 

Magnetic reconnectioii^ is a fundamental mechanism to power various events in lab- 
oratory, space, solar, and astrophysical plasmas. The reconnection process releases a fast 
plasma jet by consuming the magnetic energy stored in the upstream region. Such a plasma 
jet further interacts with the external environments and exhibits complex plasma phenom- 
ena. One of the most characteristic phenomena is a "plasmoid", a magnetic island with 
plasmas embedded inside the magnetic loop. Plasmoids are well accepted to explain the 
satellite observation in the Earth's magnetotaiP^ and morphological features in solar flare 
and coronal mass ejection^^. 

Motivated by the above and many other important reasons, the nonlinear evolution of 
magnetic reconnection and the associated plasmoid evolution has been studied over many 
decades. Significant insights have been obtained by magnetohydrodynamic (MHD) simu- 
lations in twcP^**^ and three dimensiond^^^. These simulations showed that the nonlinear 
evolution of the reconnection-plasmoid system involves complicated structures, such as slow 
shocks of the Petschek outfioW^ and various shocks around the plasmoicP^DSHSIIEl^ 

Recently, reconnection research has been rapidly developing in relativistic astrophysics, 
|-QQp3l32] |-]^jg context Zenitani et al.'^ carried out MHD simulations on relativistic mag- 
netic reconnection. They discussed several shock structures that were not known in the 
nonrelativistic plasmoid physics: the postplasmoid vertical shocks and the shock-reflection 
structure in the plasmoid- front. 

It is known that reconnection is highly influenced by upstream conditions. One key 
parameter is the Alfven speed, which usually approximates the outflow jet speed. Another 
key parameter is the plasma beta, the ratio of the plasma pressure to the magnetic pressure. 
As will be discussed later, it corresponds to the ratio of the Alfven velocity and the sound 
velocity. The beta is typically below the unity in the lower coronaP'^and the lobe (upstream 
region) of the magnetotail plasma sheetP^. 

In this paper, we explore fine structures of a plasmoid in fast reconnection when the 
upstream region consists of low-beta plasmas by using two-dimensional MHD simulations. 
Employing a high-resolution shock-capturing (HRSC) code, we successfully resolve a variety 
of shock structures. We confirm that the new shock structures in Ref. |32] similarly appear in 
the nonrelativistic regime. We further discuss the condition for these shocks, and argue that 
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they are common features of plasmoids in the low-beta plasmas. We also report some new 
features such as the Kelvin-Helmholtz-like turbulence inside the plasmoid and a potential 
signature of the corrugation instability of the slow shock surface. 

This paper is organized as follows. We describe our numerical setup in Sec. |TlJ We present 



simulation results and then investigate the local structures in depth in Sec. Ill Section IV 



presents a brief parameter survey. We discuss these results in Sec. IVl particularly focusing 



on the plasma beta condition. Section VI summarizes this paper. 



II. NUMERICAL SETUP 



We employ the following resistive MHD equations. 

dp 

di 
dpv 



+ V ■ ipv) = 0, (1) 



+V-{pvv+pTl-BB) = 0, (2) 

1^ + V ■ (^{e + pT)v - {v ■ B)B + r/j x = 0, (3) 
dB 

_ + V ■ - Si;) + V X {r]j) = 0, (4) 

where Pr = P + is the total pressure and e = p/(r — 1) + pv'^ /2 + 5^/2 is the total 
energy density. The polytropic index is set to F = 5/3. All other symbols have their 
standard meanings. 

These equations are solved by a Godunov-type MHD code. Numerical fluxes are calcu- 
lated by a Harten-Lax-van Leer (HLL) method.l^ The spatial profiles are interpolated by a 
monotonized central (MC) limiter.'^ The electric currents at the cell surfaces are calculated 
with the second order.l^S The second-order total variation diminishing (TVD) Runge-Kutta 
method is used for the time-marching.'^Sl The hyperbolic divergence cleaning method is em- 
ployed for the solenoidal condition.H^ 

The simulations are carried out in the x-z two-dimensional plane. The reconnection 
point is set to the origin (x, = (0,0). To reduce the computational cost, one quadrant 
(x > 0,2; > 0) is solved considering the symmetry of the reconnection system. We set 
a symmetric boundary at a; = 200 and a conducting wall ai z = 150. The domain of 
[0, 200] X [0, 150] is resolved by 6000 x 4500 grid cells. Harris sheet is employed as an initial 
model: B[z) = taiah{z)x, v = 0, p{z) = 1 + cosh^'^ (z) / (3up, p{z) = 0.5(/3„p + cosh~^(2;)). 
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where Pup is the plasma beta in the upstream region. Importantly, we consider that the 
plasma beta is sufficiently lower than the unity, /3„p < 1. We set f3up = 0.2 in our main 
run. By definition the Alfven velocity in the upstream region is set to CA,up = \B\/y/p = 1. 
Therefore the time is normalized by the Alfven transit time of the Harris sheet thickness. 
Note that the evolution of reconnection is best characterized by the Alfven transit time. 
The sound speed is uniform, Cs = a/TAjpT^ ~ 0.408. 

It is known that the reconnection process is critically influenced by the resistivity t]. In 
the MHD framework, a realistic fast reconnection requires a localized diffusion region in a 
vicinity of the reconnection pointP. This has been supported by simulations with spatially- 
localized resistivityJ ^"*^^ * ^^ ' Physically, such resistivity comes from local plasma instabilities. 
Hall effects,'^or electron kinetic effects,'^none of which are fully understood. It is also known 
that physics-minded anomalous resistivity models allow fast reconnection. The resistivity is 
enhanced by a phenomenological measure of microphysics, e.g., the current intensity (j/p). 
Since such quantity is the strongest around the X-point, the resistivity is localized and then 
Petschek-type fast reconnection spontaneously switches onJ ^^ * ^^ l Regarding the macroscopic 
plasmoid dynamics, it was reported that both two models basically exhibit similar results.'^ 
In this work, to keep the system simple and to see an evolution of a single reconnection, we 
employ the locahzed model, 

■r]{x, z) = + {r]i - i]o) cosh~^(Va?2~+^), (5) 

where ?7i = 1/60 is the localized resistivity and rjQ = 1/1000 is the background resistivity. A 
weak perturbation of 6 Ay = — 0.06exp[— (x^ -|- -2^)/4] is imposed to magnetic flelds to speed 
up the reconnection onset. 

III. RESULTS 

The panels in Fig. [T] present the time development of Vx, which most signiflcantly char- 
acterizes the system evolution. While we calculated the flrst quadrant {z > 0), the bottom 
half {z < 0) is inserted to understand the structures better. The evolution can also be 
seen in the movie available online. The reconnection process starts around the origin. The 
localized resistivity causes magnetic fleld dissipation there, giving rise plasma inflow from 
the upstream region. Reconnection transfers the upstream magnetic flelds toward the +x- 
direction along with an outflow jet. The flux transfer becomes faster and faster. Around 
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t > 100, the transfer rate Ey^n = rjijy^R reaches its typical value of 0.05 [cA,upBo], where the 
subscript R denotes the properties at the reconnection point. After that, the system exhibits 
a quasisteady evolution, whose reconnection speed remains fast, Ey ji ~ 0.05. Using prop- 
erties at an upstream position {x,z) = (0,5) (denoted by '), the rate can be re-normalized 
to 7^ = Ey^n/ {c'j^B') ~ 0.08. This is a fast reconnection rate of 7^ ~ (9(0.1). 

Several characteristic structures already appear in Fig. [T^. The reconnection jet becomes 
Alfvenic {vjet ~ 1) and there is a small magnetic island (plasmoid) in front of the jet. The 
plasmoid is heart-shaped^^'^, because MHD velocities are faster in the low-dense upstream 
region than ones in the plasma sheet {z ~ 0), and because dense plasmas throttle the 
outward motion of the plasmoid near the plasma sheet. One can recognize vertical structures 
of velocity jumps around a; ~ 21 and a; ~ 48 in the upstream plasmas. As the time goes 
on, the plasmoid grows bigger and bigger as it travels to the -|-x-direction. It seems that 
the macroscopic picture evolves self-similarly in later stages.'^ This is because the resistive 
spot and the Harris-sheet thickness become relatively smaller than the entire system, as the 
system grows larger. Meanwhile, fine structures become clear inside the plasmoid. Figure 
[T]d shows the profile at t = 250 of our main interest. The vertical jump structures are found 
around x ~ 60 and x ~ 100. The heart-shaped plasmoid now looks like a lobster-claw, as 
the bifurcated parts become much larger than ones in the earlier stage. 

Shown in Fig. |2^ is the plasma density p and the flow pattern at t = 250. Many plasmas 
are confined in the plasmoid. Some of the upstream plasmas are swept into the plasmoid 
by the reconnection jet. Other upstream plasmas directly enter the plasmoid across its 
outer boundary. In addition, the plasmoid eats Harris-sheet plasmas while it moves to the 
-|-x-direction. Figure [2]d presents two thermodynamic properties. The upper half shows 
the temperature including the Boltzmann constant T = p/p- In general, plasmas are hot 
inside the reconnection jet and the plasmoid. The bottom half shows a useful measure 
s = p/p^, which behaves similarly as specific entropy. Unless there is nonadiabatic heating, 
it is conserved throughout the convection d/dt[p/p^) = 0. 



In the Sees. HI A HI F , we take a closer look at local plasmoid structures at t = 250, 
from the left Petschek outflow to the front side of the plasmoid. Note that boundary effects 
are completely ruled out. Since the fastest magnetosonic speed is [c^^^p + c^]^/^ ~ 1.08, even 
the initial impact has not yet come back to the presented simulation domain. 
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A. Petschek outflow 



As presented in Figs. [T] and |2^, the fast outflow jet flows from the reconnection point. 
The jet is separated from the upstream regions by slow shocks^, and therefore one can see 
hot and high-entropy plasmas inside the outflow channel (Fig. |2}d). Figure [S] presents one- 
dimensional properties along the outflow line, z = 0. The jet speed quickly increases to 
~ 0.8-0.9, which is consistent with the Alfven speeds in the upstream sides of the Petschek 
shocks. Note that the local Alfven speeds are slightly slower than the initial speed of 1, 
because the upstream conditions evolved. The jet speed increases again in the downstream 



(x > 51), but we will discuss the downstream region later in Sec. IIIC 

To examine the shock properties in depth, we carry out Rankine-Hugoniot analysis in 
various locations in the simulation domain at t = 250. They are indicated by white marks 
in Fig. [T]d and the relevant results are presented in Table [1} We determine the shock normal 
h by the minimum variance analysid^S and check the results by considering the magnetic 
coplanarit}^. We determine the shock velocity Vsh by the mass conservation. Then, we 
compare the normal flow speed in the shock frame (i> ■ n — Vsh) with the MHD wave speeds 
in the n-direction. We calculate Mach numbers with respect to the fast-mode, the Alfven 
(intermediate) mode, and the slow-mode speeds: A^/, Aii, and Aig- The subscript 1 denotes 
the upstream properties and 2 denotes the downstream properties. The first data (No. 1) in 
Table |l] correspond to the properties across the Petschek slow shock at {x,z) ~ (40, ±1.35). 
The slope angle of the Petschek shock is small, ~0.03. It is reasonable that the angle is a 
fraction of the field line slope in the upstream side ~0.095 or the normalized reconnection 
rate TZ ~ 0.08. The shock is stationary, Vgh = 0, as expected. Across the shock, only 
the slow-mode Mach number changes supersonic to subsonic (Table |l]). This is a reasonable 
indication of the slow shock. 

Inside the outflow, we notice that there is a weak and narrow density cavity along the 
neutral line z = 0, in the closer vicinity of the reconnection point. The small box inside 
Fig. [2^ zooms up the relevant region. At x = 20, the central density is 10% lower than 
the surrounding density. This is because plasmas are hotter there, as pointed by a previous 
worliP'. The central plasmas are Joule-heated near the reconnection point, while the sur- 
rounding plasmas are heated by the Petschek slow shocks. These two may be separated by 
contact discontinuities in an idealized limilpSEZI. Since our numerical scheme (HLL scheme) 
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is weakly diffusive, the boundaries become unclear further downstream. In the real world, 
the contact discontinuities will disappear due to the field-ahgned plasma mixing. 

B. Postplasmoid slow shock 

Behind the plasmoid one can see a sharp vertical structure at x ~ 51 in Vx and T (Figs. 
|l| and (2)3). We find that this is a slow shock. We examine the shock properties near a 
typical point of {x^z) ~ (51,6). The Rankine-Hugoniot properties are presented in No. 4 
in Table |Tj To better understand the structure, one-dimensional properties along z = Q is 
also presented in Fig. |4j Note that the shock normal is virtually parallel to the x-axis. The 
normal and tangential components are similar to the x and 2;-components, respectively. In 
this case, the right side is the shock upstream, the left is the downstream, and the shock 
moves to the +x-direction at the speed of Vsh ~ +0.31. In Fig. |4^, the thin line shows the 
normal flow speed in the shock frame, f ^ = f „ — Vsh- For comparison, the local-frame slow- 
mode speed in the normal direction is presented as the dotted line. Comparing them and 
the slow-mode Mach numbers M.s in Table |T], one can see that an upstream supersonic flow 
in the right side becomes subsonic in the left downstream. The shock exhibits many other 
slow-shock features. Along with the compression (Fig. |4|d) and the adiabatic heating (T in 
Fig. |2}d), one can see that the tangential magnetic field Bt (~ Bz + 0.03) decreases across 
the shock (Fig. |4|d). Meanwhile, this is a nearly parallel shock and magnetic dissipation 
is quite limited. The measure s is virtually unchanged (Fig. ^p) and the heating is weak 
{T2/T1 in Table |ll). 

As discussed in Ref. |32l this vertical shock is related to a reverse plasma flow around 
the plasmoid. When the plasmoid moves to the right, it compresses plasmas in the nearby 
upstream region. Consequently, gas pressure drives plasma flows in the field-aligned direc- 
tions. In addition, since there is a density cavity behind the plasmoid (Fig. [2^), the reverse 
flow in the postplasmoid region (50 < x < 70) is particularly strong, as can be seen in Fig. 
[T]d < 0; blue regions). The reconnection inflows usually move to the ±2-directions, but 
they are accordingly modulated by the reverse flows (Fig. [2^). Due to the field-aligned 
expansion, the temperature becomes lower in the reverse flow region (Fig. ^p)- The slow 
shock is located at the front, where the reverse flow supersonically hits the right-going shock 
front. 
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C. Postplasmoid outflow 

We find that the postplasmoid vertical shocks affect the reconnection jet inside the outflow 
channel. The vertical shocks meet Petschek shocks at x ~ 51.2. In Fig. [3} one can see that 
the reconnection jet suddenly becomes faster around there. The jet further travels in the 
+x-direction, and then it hits the plasmoid at a shock at x ~ 61.2. The results of the 
shock analysis are shown in No. 3 in Table |l| This is a reverse- type fast shockPI^l. The 
shadow in Fig. [3] indicates the outflow region between the vertical shock and the fast 
shock (51.2 < X < 61.2). It seems that the pressure gradient accelerates the jet in this 
shadow region. The jet speed reaches = 1.2 at the fast shock, which is much faster 
than the initial upstream Alfven speed, CA^up = 1. A similar effect is reported by Shimizu & 
Ugg^jHUig]^ who claimed that deconfinement of the outflow channel leads to an adiabatic-type 
fluid acceleration. 

We think that two mechanisms work in the shadow region. First, the system adjusts 
the balance condition of the Petschek shock, based on the upstream conditions. In the left 
quiet region {x < 51.2), plasmas from the downstream side of the vertical shock come in 
to the Petschek shock. The outflow speed is well approximated by the upstream Alfven 
velocity there. In the shadow region (51.2 < x), plasmas from the upstream side of the 
vertical shock come in to the Petschek shock. The upstream Alfven speeds are faster there, 
CA,up' ~ 1.16, mainly due to the lower density (Fig. |4|d). Then, the Lorentz force at the 
Petschek shock accelerates plasmas to the "new" outflow speed of ~CA,up'. One can see a 
rarefaction- like transition between the left Petschek flow and the right new Petschek flow in 
Fig. [3j Second, an adiabatic-type acceleration further boosts the jet.'^^'^ This is a simple 
hydrodynamic effect. Since the outflow speed exceeds the local sound speed of ~ 0.6-0.7, 
such a supersonic flow is adiabatically accelerated when its cross section increases like a 
nozzle. It particularly works near the fast shock, where the plasmoid opens up the outflow 
channel and so the final outflow speed of ~ 1.2 exceeds the maximum upstream Alfven 
speed of CA,up' ~ 1-16. 

According to the Rankine-Hugoniot analysis (No. 2 in Table [l]), the Petschek shock in 
the postplasmoid region is stronger than the left-side (No. 1). The shocked plasmas are 
hotter and the entropy is higher than plasmas which come from the left Petschek shock, 
as seen in Fig. ^jp. A sharp boundary separates plasmas from the left and right Petschek 
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shocks. We think this boundary is a contact discontinuity, which originates from the shock- 
shock interaction point. Similarly, we expect that the field-aligned mixing will blur the 
discontinuity in the real world. 

D. Outer shell and external slow shocks 

The plasmoid is surrounded by a relatively sharp boundary. The Rankine-Hugoniot 
analysis across this outer shell (Nos. 5 and 6 in Table [l| tells us that it is a slow shockf^. 
Judging from the heating ratio of T2/T1, the shock is stronger in the back side of the 
plasmoid than that in the front side, as discussed in previous literature— . In addition to the 
reconnection jet, the slow shock on the outer shell supplies hot plasmas into the plasmoid. 

Outside the plasmoid, we find another pair of slow shocks around {x,z) ~ (100, ±10). 
The Rankine-Hugoniot result is presented in No. 7 in Table |Tj This shock is rather weaker 
than the other shocks. In fact, one can hardly find a jump in the measure s in Fig. [2)3, 
because dissipation is quite limited. However, our code does capture a clear jump structure 
in the velocity (Fig. [T|3). Interestingly, as can be seen in Fig. [T^ and the movie available 
online, the shock is originally located near the front edge. It gradually shifts backward to a 
stable location as the system evolves. 



E. Internal structure 



The reconnection jet hits the plasmoid at the fast shock (Sec. HI C), where plasmas enter 
the plasmoid. In the downstream side, there is a busy transition region at x ~ 70, where the 
shocked fiow further slows down and accumulate the magnetic fields Bz (Fig. [s]). This region 
corresponds to the loop-top fronlP^l. Plasmas hit the loop-top front and then they diverge 
from the central region to the upper and lower regions inside the plasmoid. The reconnected 
fields are piled up between the loop-top front and the sharp tangential discontinuity (TD) 
at X ~ 80. The field is remarkably strong in the left side of the TD. The TD separates 
upstream-origin plasmas and plasmas initially located inside the Harris current sheet. This 
is visible in the s-profile in Fig. ^jp. The entropy is high in upstream-origin plasmas, because 
they are shock-heated or Joule-heated. On the other hand, the right plasmas just come from 
the low-entropy Harris sheet. 
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Figure [5^ presents the density profile inside the plasmoid. While the plasmoid is dynam- 
ically growing and expanding, the TD travels in the +x-direction at a nearly constant speed 
of Vt ~ +0.5. To understand the plasma motion better, the velocity profile in the moving 
frame of Vt is shown by arrows. Plasmas at 2; ~ constantly travels in the — a;-direction 
from the right. After they hit the TD, they turn to the +z direction, and then they are 
reflected to the +x-direction around ^ ~ 2. The TD works as a "magnetic wall" to reflect 
plasmas from the right side. In the original frame, this traveling wall hits the immobile 
Harris-sheet plasmas (f ~ 0). Therefore, the reflected flows are twice faster than the wall 
motion in Fig. [TJo, ~ 2vt ~ 1. 

The fast reflected flow triggers velocity-shear-driven turbulence. In the online movie, one 
can see that a vortex rolls up where the flow turns its direction and that the reflected flow is 
deformed to a kink-type structure. We think this is due to the Kelvin-Helmholtz instability 
between the reflected jet and the surrounding medium. For example, typical parameters 
on the upper side of the reflected jet are B ~ Bq/2 ~ 0.5, p ~ p(t = 0,z = 0) ~ 6, and 
Av Vt ^ 0.5, where Av is the velocity shear. From the linear theory,'^^! we find that 
the shear is strong enough to bend the magnetic field lines, Av > 2\B\/^, but is slower 
than the sonic upper limit, Av < 2cs. Therefore, the shear layer would be unstable to the 
Kelvin-Helmholtz instability. One can see that the reflected flow is violently flapping in Fig. 
[T]d and the online movie. The reflected flow and the resulting turbulence mix plasmas inside 
the plasmoid. 

In Fig. [5^ and the online movie, one can also see that the incoming flow repeatedly flaps 
in the right region of 98 < x and z ~ 0. This is related to the shock-reflection structure 
which will be discussed in Sec. IIII Fl 

F. Oblique-shock reflection 

Shown in Fig. |5|d is the fa.-profile in the front side of the plasmoid. The outer shells 
reach the leading edges of the plasmoid at (x, z) ~ (121.8, ±2.4). Around the central region, 
there is immobile plasmas (f ~ 0), because high-density Harris-sheet plasmas are initially 
confined here. Outside the center, plasmas fast travel in the +x direction. These twin flows 
are essentially driven by the reconnection outflow or the reflection by the magnetic wall. The 
central plasmas and the twin flows are separated by transition regions. The magnetic fields 
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and plasmas flows are rather parallel to the transition regions, but they gradually enters 
from the central Harris sheet to the plasmoid flows. 

As pointed out by Abe & Hoshino'^, the transition region looks like an intermediate 
shock, across which the tangential magnetic flelds change their polarities. Unfortunately, our 
Rankine-Hugoniot analysis fails to tell whether or not this is an intermediate shock (No. 8 
in Table [l]). This is probably due to the following three reasons. First, the assumption for the 
analysis may be invalid here. The Rankine-Hugoniot analysis deals with a time- stationary 
and one-dimensional structure. However, oblique shocks repeatedly hit the boundary regions 
as discussed later. A time-dependent flow and a weak curvature makes our analysis difficult. 
Second, the slow and intermediate speeds are hard to separate when the field lines are 
quasiperpendicular to the normal direction. Third, our numerical scheme (HLL scheme) 
is not ideal for the structure here. The HLL code excellently captures the fast shocks but 
handles other discontinuities in an upwind manner. When the characteristic speeds are slower 
than the fast-mode speeds by order of magnitude, the code smoothes the structure. Here, 
the fast speeds are typically ~ 0.5-0.6, while intermediate and slow speeds are on an order 
of ~ (9(0.01) [cA,Mp]. These three reasons make it difficult to discuss the slow/intermediate 
Mach numbers. In Table [T], we add asterisk signs to such unreliable Mach numbers. The third 
issue will be improved by a better scheme which fully considers intermediate discontinuities. 
In the right side near the flow edges (No. 9 in Table |l]), the boundaries look like slow 
shock^. Unfortunately, we fail to classify this shock for the same reasons. We expect that 
the intermediate shock changes to the slow shock through the switch-off shoclP^'^. 

In the central channel, one can see that the oblique shocks propagate from the twin flow 
edges (Fig. |5}d). They are fast shocks (No. 10 in Table |l]). This is because the edges travel 
faster than the fastest right-going wave inside the central channel. These oblique shocks 
are repeatedly reflected by the intermediate shock region^. Such shock-reflection leads to 
a chain of the diamond-shaped patterns, as denoted by the "diamond-chain" structure in 
Ref. |32l The second and third right ones look irregular in this case, but one can see that all 
other rectangles are excellently in order in Fig. [5]d. The repeated shock-reflection explains 
a violent fluctuation in sonic properties in the front region (95 < a; < 120) in Fig. |3] 

Roughly speaking, the plasmoid edges travel at a sizable fraction of the reconnection jet 
speed Vjet, which is well approximated by the initial upstream Alfven speed, CA,up- On the 
other hand, since magnetic flelds are very weak inside the current sheet, the right-going 
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fast-mode speed is approximated by the sound speed in the current sheet. Therefore, a 
requirement for the obhque shock and the diamond-chain is 

Cs,cs ^ CjX^upi (6) 

where the subscript cs denotes the current sheet properties. Needless to say, our initial 
configuration satisfies this: Cg^cs = Cg = 0.408 and CA,up = 1- 

In the fronst side, there are small bow-shock-like structures around x = 121.8, z = 2.55- 
2.75, extending outward from the plasmoid edges. We do not recognize steep shocks there. 
We think they are the wave fronts of the fast modd^^l. 



IV. BETA DEPENDENCE 

We carry out another two runs using different f3up parameters, f3up = 0.1 and 0.5. All 
other parameters are the same as the main run. Figure |6^ represents a snapshot with a 
lower upstream beta, /3„p = 0.1. The system builds up slower than the main run, because 
the system has relatively denser plasmas in the downstream Harris sheet. One can recog- 
nize many characteristic signatures except for the forward vertical slow shocks outside the 
plasmoid. Actually, there are slow shocks near {x,z) ~ (82, ±9.8), but they are too small 
to recognize in the figure. The postplasmoid vertical shocks are located at a; ~ 40. The 
velocity jump at the shock is smaller than the main run, because the sound speed is initially 
slower, Cg = \/YP>upl'2. ~ 0.29. 

Figure [6|d shows the other run with f^^p = 0.5. One can see that the postplasmoid slow 
shocks are located in the back side of the plasmoid's outer shell, a; ~ 55. The other pair of 
slow shocks are located near the front edges of the plasmoid, x ~ 115. At a glance the front 
shocks look like switch-on fast shocks, but they are switch-off type slow shocks. Magnetic 
fields are gradually bent in the right side of the forward vertical slow shocks. Interestingly, 
a finger-like structure can be seen on the slow shock surface, around 2 < 1^1 < 6 outside the 
plasmoid- front. It slowly grows there and it may eventually destroy and blur the forward 
slow shocks. We think that this is a signature of the corrugation instability of the slow shock 
surfac^. However, detail analysis of the instability is beyond the scope of this paper. 
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V. DISCUSSION 

In this paper we visited many structures in and around the plasmoid. They are schemati- 
cally illustrated in Fig. [7j This extends an informative result in an earlier work (i.e., Fig. 10 



in Ref. ITSll . Among these structures, the postplasmoid vertical shocks (Sec. IIIB) and the 



shock-reflected diamond-chain structure in the front side (Sec. IIIF) were recently found in 
the relativistic regime.l^Sl We conflrmed them in the nonrelativistic regime as well — they are 
ubiquitous, regardless of the relativistic effects. A reverse flow or a backward flow around 
the plasmoid is often found in recent large-scale simulation^^^'^, however, we clearly demon- 
strated that a slow shock stands at the front of the reverse flow. We also introduced another 
vertical slow shock outside the plasmoid. Those two vertical shocks are similarly found in 
other runs in our parameter survey. 

We think that the vertical slow shocks are logical consequences of the low-beta condition 
of f3up < 1. From the following relation, 

ics^up/cA,up) = (r/2)/3„p, (7) 

we see Cg^up < CA,up when < 1. As already discussed, the plasmoid-front propagates at 
a sizable fraction of the reconnection jet speed vjet ~ CA,up- When f3up < 1, the plasmoid 
supersonically travels in the system. At some location, the expanding plasmoid structure 
interacts with the surrounding plasmas at a supersonic speed, and then a shock stands 
there. Another interpretation is as follows. In a self-similarly expanding coordinate of the 
plasmoid system,'I3the surrounding plasmas supersonically flow in the — a;-direction near the 
plasmoid-front, while the relative flow is very small near the reconnection point. A shock 
stands at the point where the supersonic flow slows down to subsonic speed. In summary, 
when Pup ^ 1, the plasmoid system will inevitably involve vertical slow shocks. 

In addition, the compression effects by the plasmoid introduces another pair of vertical 



shocks outside the plasmoid. As discussed in Sec. IIIB the compression invokes the fleld- 
aligned flows in the ±x directions. In addition, the adiabatic heating or cooling weakly 
modify the local sound speed, Cg = VTs^" ■ oc p^^^. In the main run, the relative 

velocity between the plasmoid system and the upstream plasmas are initially supersonic 
in the right side. It temporally becomes subsonic at the forward vertical shock x ~ 100 



(Sec. HID) due to the compression effects. However, due to the fleld-aligned acceleration in 
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the — x-direction, the flow becomes supersonic again, and then it becomes subsonic at the 



postplasmoid vertical shock at x ~ 51 (Sec. IIIB). 

The positions of the vertical shocks depend on the Pup parameter, because it controls 
the ratio of the sound speed to the Alfvenic plasmoid motion. When Pup is lower than the 
main run {Pup = 0.1), the expanding speed of the entire structure easily exceeds the sound 
speed. Therefore, the postplasmoid vertical shock is closer to the reconnection point, as 



demonstrated in Sec. IV (Fig. [6fi). The forward vertical slow shock can be small, when 
the plasmoid compression effects are not strong enough to change the supersonic condition. 
Probably the forward vertical shock disappears in the extreme limit of Pup ^ 1. On the 
other hand, when Pup is higher [Pup = 0.5), the vertical slow shocks are nearer to the 
front-side of the plasmoid: the postplasmoid and forward ones are located just behind 
and in front of the plasmoid, respectively. As Pup increases, the postplasmoid slow shock 
moves frontward even in the backward side of the plasmoid outer shell, and then it will 
eventually disappear. The forward slow shock will disappear when the Alven speed is no 
longer supersonic. Pup ~ 2/T = 1.2. Depending on /3„p(< 1), we expect one or two pairs of 
vertical slow shocks around the plasmoid. 

Our discoveries can readily be applied to various numerical results in low-beta plasmasJ ^^ * ^° * ^^ * ^^ ^ 
Abe & Hoshino'^ carried out large-scale MHD simulation of a plasmoid in a similar con- 
figuration. While they did not mention it in the text, one can see a shock-like structure 
outside the plasmoid in RUN 2 with Pup = 0.25 (the top panel in Fig. 9 in Ref. [18]). We 
think that this is a forward vertical slow shock. Zenitani et al.^ recognized a small shock 
outside the plasmoid (indicated by an arrow in Fig. 2c in Ref. [321) in the relativistic regime 
with Pup = 0.1. Inspired by this work, we carry out further analysis and confirm that it was 
a forward vertical slow shock. In multiple plasmoid systems, Tanuma & Shibata'^ found 
"bow shocks" during the nonsteady reconnection under an anomalous-type resistivity with 
Pup = 0.2 (Fig. Id in Ref. [201). Judging from displacements, their plasmoid velocities are 
sub-Alfvenic, on an order of ~Cs,«p < ca,up- In such a developing stage, the forward vertical 
shock is located near the front, and so it would be a forward slow shock. Finally, a recent 
work by Yu et al. reported a steepening structure of slow-mode wave behind the primary 
plasmoid (Fig. 15a in Ref. [2^. It would be related to a postplasmoid slow shock, as we 



found it in our relevant run with Pup = 0.5 (Sec. IV Fig. |6|; 

The shock-refiection or the diamond-chain structure is also a characteristic in the low-beta 
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condition. Dropping a factor of order unity, the requirement (Eq. Ig]) reads 



^s,cs 



CA 



%p<l. (8) 



Although we still have a free parameter {cs^cs/cs,up), the condition can be easily satisfied in 
the \ow-f3up regimes. 

Let us discuss some limitations of our simulation model. First, we assumed that the 
evolution in the bottom quadrant {z < 0) is symmetric to that in the upper quadrant 
{z > 0). We do not think the quadrant assumption affects the main shock structures, 
because Ref. [32] was carried out without the assumption. On the other hand, we found 



the Kelvin-Helmholtz-like turbulence in the internal region (Sec. HIE). Since this class of 
instabilities may prefer odd-parity perturbation across the neural line {z = 0), the central 
Harris sheet may flap in very late stages. This is left for future investigation. Second, in 
the three dimensions, it is known that quasiperpendicular slow shocks are unstable to the 
corrugation instability.'^ We do see a potential signature of the instability even in the two 



dimensions (Sec. IV). In the real world, the vertical slow shocks may be violently folded and 
then we will see turbulent transition layers instead. Regarding the shock-reflected diamond- 
chain structure, the oblique fast shocks will be stabl^, but it is not clear how stable the 
intermediate shocks are. It is interesting to see how those structures are modulated in the 
three dimensions in future. 



VI. SUMMARY 

We carried out HRSC MHD simulations to study the fine structures of a plasmoid. We 
introduced new shock structures: (1) postplasmoid vertical slow shocks, (2) forward vertical 
slow shocks, and (3) shock-reflection structures, some of which were recently reported in Ref. 
[321 We further found new features such as (4) the flapping motion of the reflected jet and 
(5) the finger-like instability of the slow shock surface. We interpreted that the vertical slow 
shocks and the shock-reflection are consequences of the Alfvenic plasmoid motion, which 
is supersonic in low-beta plasmas. We argue that these shocks are rigorous features of a 
reconnection system in low-beta plasmas. 
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(b) Temperature and entropy 10^ 




FIG. 2. (Color online) Snapshots at t = 250, overlayed by magnetic field lines, (a) Plasma density 
p. Arrows indicate plasma velocity: \v\ = 1 is projected to the length of 5. The inner box shows 
the smaller doamin of x G [15, 25] and z £ [—1.5, 1.5]. (b) Plasma temperature T in the upper half 
and the entropy measure s = p/p^ in the bottom half. 
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FIG. 3. One dimensional cut of physical properties along the neutral line, z = 0. The reconnected 
magnetic field B^, outflow velocity v^, plasma pressure p, and the density p. The shadow indicates 
the region between the vertical slow shock {x = 51.2) and the fast shock {x = 61.2). 
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FIG. 4. (a) One dimensional cut along z = +6. Normal velocity Vn (~ Vx), tangential velocity 
Vt, (~ Vz), shock- frame normal velocity v'^, and the shock- frame slow-mode velocity in the —x- 
direction. (b) Tangential magnetic field Bt, pressure p, and density p. See the second vertical axis 
to measure p. 
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FIG. 5. (Color online) (a) Plasma density p inside the plasmoid at t = 250. The white lines 
indicate magnetic field lines. The arrows indicate the plasma velocity in the moving frame of 
Vx = +0.5. The spatial length 2 corresponds to \v\ = 1. (b) Profile of Vx at t = 250 in the front 
side. (Enhanced online) 
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(a) P - 0.1 275.0 xio° 




FIG. 6. (Color online) The f^-profiles in different simulation runs: (a) = 0.1, t = 275 and (b) 
Pup = 0.5, t = 210. 
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. / 10 




FIG. 7. Schematic drawing of plasmoid structures: 1. Petschek slow shoclP, 2. outer shell (slow 
shock)ff^, 3. intermediate shocls?^ and slow shock^^, 4. fast shodi^^'^, 5. loop-top fronlP^, 6. 
tangential discontinuity, 7. postplasmoid vertical slow shoclP^l^ g. forward vertical slow shock (Sec. 



HID), 9. fast-mode wave front^^ , 10. shock-reflection (diamond-chain)'^-^, 11. contact discontinuity 



(Sec. IIIC), A. reconnection inflow, B. outflow jet, C. postplasmoid backward flow, D. internal 



flow, and E. flapping jet. 
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